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Abstract 

We estimate the cross section for 7p — > Tp by two independent methods. First, by 
studying the corrections to the naive leading-order QCD formula and, second, by using 
parton-hadron duality. The estimates are in good agreement with each other and with 
the recent measurements of the cross section at HERA. 



The elastic photoproduction of T mesons has recently been measured for the first time 
both by the ZEUS Jl| and HI collaborations. The observed cross section Bo = <t(tp — > 
Tp) B(T — > is found to be rather large 

Ba(T) = 13.3 ± 6.0 l^pb at (W) = 120 GeV (ZEUS [1]), (1) 



Ba(T) = 16.0 ± 7.5 ± 4.0 pb at (W) = 160 GeV (HI (preliminary) @) (2) 

where T stands for the sum over the T(IS'), T(2S") and T(3S) states and where (W) is the 
average photon-proton centre-of-mass energy. In both presentations it was noted that the 
measurements are higher than the predictions of various QCD models by about two standard 
deviations. 

In this paper we estimate elastic T photoproduction from perturbative QCD by two inde- 
pendent methods. Both calculations are found to give cross sections which agree well with the 
data, and so we conclude that there is no disagreement with QCD. The first method is to study 
the main corrections to the naive QCD formula for heavy vector meson photoproduction. 
The second method makes use of parton-hadron duality similar to where p electroproduc- 
tion at HERA was described successfully. 



Method I 

The first method is based on the leading-order expression for the photoproduction of a heavy 
vector meson of mass My HI 
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where Y v ee is the partial width of the V — > ee decay, as is the QCD coupling, a = 1/137 is 
the QED coupling and g(x,fi 2 ) is the gluon density measured at x = M v /W 2 and the scale 
H = ttiq ~ My/2. We discuss below four types of corrections to this naive formula. 



(a) Relativistic corrections 

A non- relativistic wave function for the heavy vector meson was used to derive (|^), so first 
we quantify the size of the relativistic corrections. For J/ip photoproduction these corrections 
have been the subject of debate || ||, [7]]. Since the velocity of the charm quarks can be sizeable 
in the J /if) meson ((f 2 ) — j\ the kinematical correction coming from the quark propagators 
suppresses the J /if) cross section by a factor of about 2 according to Ref. or up to a factor of 
10 in [7j. However it was shown in that to a good approximation this correction is nullified in 
the cross section formula (Q) since it has been written in terms of the J / if) mass, instead of the 
current quark mass (2m c < which should have been used in perturbative QCD. So it turns 
out that there are practically no (v 2 ) corrections to the leading order result (|3|). This result is 
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consistent with a more detailed study by Hoodbhoy ||. Strictly speaking to estimate the (v 2 ) 
corrections we should also take into account the contributions coming from more complicated 
components of the vector meson wave function which contain one or more gluons besides the qq 
pair. These components were included in a self-consistent way by Hoodbhoy M, who showed 
that the total 0(v 2 ) relativistic correction (including corrections to the quark propagator and 
extra gluons in the J j ip wave function) to the naive cross section formula (||) amount to at most 
only 7% for J ftp photoproduction. The correction should be even smaller for T production. 
Therefore we neglect these relativistic effects below. 



(b) Real part 

Only the imaginary part of the '-/p — > Vp amplitude is expressed in terms of xg(x,p 2 ) in 
(|3|). To restore the real part we may use the approximation 

Re A/lm A ~ vrA/2 (4) 

where 

A = d log A = 1 d(xg(x,p 2 )) ^ 
Slogs xg(x, p 2 ) dln(l/x) 

This approximation is valid for small A, say A < ~, for the even signature amplitude A. The 
real part correction factor 

/ 7T 2 A 2 \ 

C b = (1 + -^-J = 1.43(1.36) or 1.54(1.46) (6) 

for (W) = 120(160) GeV and scale p 2 = 25 GeV 2 , according to whether the MRS(R2) @ or 
the GRV |1(J gluon distribution is used. 



(c) The effect of off-diagonal partons 

At leading order the elastic photoproduction of vector mesons is mediated by two-gluon 
exchange. Strictly speaking the amplitude is proportional to the off- diagonal gluon distribution 
x'g(x, x'), where the momentum fraction x' carried by the "second" gluon is much smaller than 
the fraction x ~ M v /W 2 carried by the "first" gluon. The effect, R = x'g(x,x')/xg(x), was 
estimated in Ref. The correction is almost negligible for p production, giving about 10% 
enhancement of the J/ip amplitude, but is much more important for T photoproduction since 
the scale p 2 ~ M\ and the value of x = M\jW 2 are much larger. In fact the correction factor 
to the cross section formula (0) for T production is [IT] 

C c = R 2 ~ (1.4) 2 ~ 2. (7) 



(d) NLO corrections 

Finally we have NLO corrections which come from an explicit integration over the loop 
which corresponds to the convolution of the exchanged gluons with the vector meson wave 
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function. The integral is not of a pure logarithmic form, which was assumed in deriving the 
leading order result (||). The typical range of integration corresponds to the scale /i 2 = \M V and 
the non-logarithmic contribution is well approximated by a NLO correction with a coefficient 
of about | |J. That is the correction can be written as 

C d ~ (l + a s (\M 2 ) /2) ~ 1.2 (8) 

for T photoproduction. 

We are now in the position to estimate the cross section for elastic T photoproduction. If we 
assume that the elastic T differential cross section, da/dt ~ exp(frt), has a slope 5 = 4 GeV -2 , 
as for elastic J / ip photoproduction, then the leading order formula gives 

Ba(T ls )\ LQ ~ 1.7 pb (9) 

at (W) = 120 GeV. Taking into account the corrections Cb, c ,d we obtain 

Ba(T ls ) ~ 6pb. (10) 

However the cross sections measured by the ZEUS and HI collaborations, ([[]) and @, did 
not select the pure T(15) state, but rather the data were integrated over an interval of the 
p, + p~ mass which includes at least the IS, 2S and 35 resonances. We estimate the combined 
contributions of the 2S and 3S states to be 40% of the T(1S) cross section, in agreement with 
the measurement of the CDF collaboration [|12|]. Thus multiplying ([TDD D Y a factor of 1.4 we 



estimate 

Ba(T) = 8.4 pb, (11) 

which is compatible with the ZEUS measurement given in (H). Going from (W) = 120 GeV to 
160 GeV, the cross section is predicted to increase by a factor of 1.6, due to the larger gluon 
density at smaller x, and so the prediction is in agreement with the preliminary HI measure- 
ment quoted in (0). 



Method II 

The second method that we use to estimate the elastic photoproduction of vector mesons 
from perturbative QCD is based on parton-hadron duality. The procedure [|J is to calculate 
the amplitude for open qq production, then to project the amplitude onto the J p = 1~ qq 
state and finally to integrate the cross section over an appropriate interval AM of the mass 
of the qq pair which includes the resonance peak. As there are almost no other possibilities 
for hadronization at M q q ~ My we should obtain a reasonable estimate of the cross section 
for vector meson production. This framework was found || to describe successfully the energy 
and Q 2 dependence of p- meson electroproduction (7*p — > pp), both for longitudinally and 
transversely polarized p mesons, including the Q 2 dependence of the <Jl/&t ratio. 

The computer code used here to estimate the 7p — ► Tp cross section takes into account 
all the corrections mentioned above, with the exception of the enhancement due to the use of 
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off-diagonal gluons. (The enhancement is negligible for p production.) We therefore include the 
extra factor C c ~ 2 of (0). As the recent measurements ([!]) and (0) cover a rather large p + pT 
mass interval^ (AM ~ 2 GeV) over the T resonances a description based on parton-hadron 
duality might even be more appropriate to describe the data. 

In practice we have to extend our code to include the mass m of the heavy b quark. First 
we restore the mass term mSw in the expression for the 7 — > qq matrix element, see Eq. (32) 
of Ref. [pl|] . Here we will use the notation A = i = +, — and A' = %' = +, — for the helicities of 



the quark and antiquark. The cross section is then given by (see also 

d 2 a T n 2 e 2 a 

t=0 = 3(Q 2 + M 2 ) 



dM 2 dt 



J 2dz \ B ii'\ 2 (12) 



where the helicity amplitudes are (for photon helicity +1) 



B ++ = B — = (13) 

2jz(l-z) 



and 

= ~ ZkTlT , = (l-*)*^. (14) 

'z(l-z) Jz(l-z) 



The variable z is the momentum fraction carried by the quark and kx is its transverse momen- 



tum. The integrals II, It are defined in [13]. The second modification due to the mass of the 
quark is associated with the fact that the 7 — > qq helicity amplitudes B^i are defined in the 
proton rest frame (pRF), while the projection onto the J p = 1~ state was done using the quark 
helicities in the qq rest frame (qqKF). Unfortunately helicity is not a good quantum number 
for a heavy quark. It can be changed by a Lorentz boost. So we have to compute the helicity 
amplitudes Ajf in the qqKF in terms of B^i in the pRF. We have 

A/i' = X! Ci i Bu' (15) 

where we will calculate the coefficients (c^ for the quark and Cj^i for the antiquark) via the 
polarized quark density matrix p. For a quark with 4-momentum and polarisation vector a M 
(satisfying a 2 = —1 and a.k = 0) we have |H 



p = Oi + m) M (16) 

where (1 + 75^/2 projects onto the state with polarisation vector a M . For the states with 
helicities j = ± in the qqKF these vectors take the form 

= (a ;a T ,a 2 ) = (a ; a) = ±(k;k k/k)/m (17) 



1 Note that in the parton-hadron duality approach AM is related to the formation time rather than to the 
width of the resonances. The 775 — > Tp data are integrated over an interval larger than the formation time but 
which includes at least the first three T states and possibly some small contribution from the bb continuum. 
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where k = \k\. Let 6^ be the analogous polarisation vectors describing helicities i = ± in the 
pRF. After the boost from the pRF to the qqKF these vectors are given by 

6* - ± (— - —■ - (z- ^ + — ^ risi 

" ~ ^2m zAf' m' m V 2 2 J + zM j W 

where M is the mass of the qq pair, z = | + kcos9/M and # is the quark decay angle in the 
ggRF. 

Now that we have all the polarisation vectors in the ggRF we can calculate the product of 
the two different projectors^ 



iTr[(l + 75^)(l + 7 5 r)] = l-(o.6) = 24. (19) 



Thus we have 



1 - (a.b) 



C++ = c— = c+_ = -c_+ = d ^— , (20) 



where c h is chosen to be negative to ensure the orthogonality of the h l states when expressed 

in the a? basis. Note that if the quark mass m — > then 

aj - &J = ± (21) 

leading to the unit matrix = 5^. 

The procedure is therefore to decompose the original amplitude Ba> in the pRF in terms of 
the amplitudes A,y with helicities = +, — in the qqRF, as in (|I~5l). The resulting helicity 
amplitudes are projected onto the J p = 1~ qq states so as to obtain the production amplitudes 
in spin states |1, m) with m = 0, ±1 

T ™ = E Vl /! dcose^vdlje)^,^). (22) 

The cross section of diffractive J p = 1~ qq pair production is then obtained by summing over 
the m = 0, ±1 amplitudes squared in an analogous way to that described in Ref. |Q. 

The expressions for 1^ of ( Jl3| ) and ( Jl^ ) were written in for the (dominant) imaginary 
part of the amplitude. They are integrals over the gluon transverse momentum £t with weight 
w corresponding to the gg-loop 

JmI L , T = J d£ 2 T f(x,e 2 T )w L>T (4,Q\M 2 ---) (23) 

where / is the unintegrated gluon distribution 



= dx9 tt\ T> - (24) 



We omit the projector A($ + to) of (M) so that the trace includes both quark and antiquark components. 



Hence the factor 2 in 2cfj on the right hand side of ( p_9| ) 



The real part of the amplitude has been included using the same approximation as in Eqs. (|4], 
U), by replacing f(x, with the derivative (tt /2)df(x, £^)/d\n(l/x), that is we have computed 

Re/L,T = \ J dl% w L>T (f T , Q 2 , M 2 , ...). (25) 

Finally the NLO corrections (analogous to Cd of (|8])) are approximated in the code by use of 
a X-factor, see [EL O. As in |13| (and (B)) the scale in the coupling as(/z 2 ) was chosen to be 



fi 2 = M 2 /4. 

To compare with the ZEUS data for 7p — > Tp we integrate over the mass interval 8.9- 
10.9 GeV. We take the mass of the b quark to be mj, = 4.6 GeV and, as in Ref. j|, we use the 
MRS(R2) set of partons ||. In this way we find 

°(lP^Y. T P) - 620 P b ( 26 ) 

at (W) = 120 GeV. This should be compared with the value 635 ± 310 pb obtained from the 
ZEUS measurement. The "ZEUS" value is calculated from (|l|) by assuming the ratios 

Ba(T ls ) : Ba{T 2S ) : Ba(T 3S ) ~ 0.7 : 0.15 : 0.15, (27) 



which are consistent with the expectations of (|3|) and with the CDF data [121, and using 



the known leptonic branching ratios of the T states. What is the uncertainty in our cross 



section estimate ( p6|) due to the choice of m^? The PDG [|15[ give 4.1 < m& < 4.4 GeV for 
the running mass evaluated at \i — nib in the MS scheme, which corresponds to the range 
4.5 < rrib < 4.8 GeV for the pole mass. The latter mass is relevant for our perturbative 
calculations. If we were to use = 4.4 or 4.8 GeV then the cross section ( p6|) would become 
900 or 380 pb respectively. 

Fig. 1 shows the ZEUS and HI data, together with our QCD estimates obtained from the two 
independent methods. The comparison is given for the 7p — > T(lS)p cross section as a function 
of the photon-proton centre-of-mass energy W. To obtain the T(IS') cross section for method 
II, in which we integrate over the mass interval 8.9-10.9 GeV, we divide the 7p — > (J2 T)p 
result by a factor 1.7. This factor is obtained using For completeness we also show in 

Fig. 1 the prediction from the naive leading order formula of method I as well as the result 
of method II without the corrections from the real part of the amplitudes. We note that the 
7P — > Tp process is particularly sensitive to off-diagonal gluon effects. These effects give an 
enhancement of about a factor of 2, which appears to be required by the data. 

To estimate the ratio of T to J/ip photoproduction we calculate both cross sections using 
all the corrections listed for method I. For (W) = 120 GeV we find 

<7(Tis)/<7(J/V) ^ 1/400 (28) 

where the MRS(R2) gluon has been used. This is near the lower bound observed by ZEUS 
[|1| . To be specific we underestimate T production and overestimate^ J / if) production. We note 

3 Note that in Ref. || absorptive corrections were included and that a better description of J/ip photoproduc- 
tion was obtained. Here absorptive corrections were not discussed as they should be small for T photoproduction. 
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that if, for example, the steeper GRV gluon was used the ratio would be 1/1500. Clearly, as 
expected, the ratio is very sensitive to the x and scale dependence of the gluon in the region of 
very small x and low scale that is sampled by J / ip production. In the case of T production we 
sample x ~ 0.01 and Q 2 ~ 25 GeV 2 where the predictions are much more stable with respect 
to different parametrizations of the gluon. 

We conclude that the recent measurements of elastic T photoproduction are in good agree- 
ment with the expectations of perturbative QCD. Of course when more precise data become 
available the effects that are estimated here — relativistic corrections, NLO contributions, off- 
diagonal parton effects and the contribution of the real part of the amplitude — should be 
calculated in detail. The elastic photoproduction of J/ip and T at HERA will be a good labo- 
ratory to probe these effects. 
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Note added 

While writing this note we received a paper on the same subject by Frankfurt, McDermott 



and Strikman fll6fl . Their procedure is comparable to our method I and also yields a cross 
section in agreement with the HERA data. However although our final values are similar, there 
are large, physically significant, differences in the various component factors. Unlike the present 
work, they do not make use of the result of Hoodbhoy M and so have a large suppression from 



the relativistic corrections. In [16] this suppression is compensated by a larger off-diagonal 
effect (2.6 as compared to 2), a larger real part and an enhancement due to re-scaling (giving 
an effective Q 2 of about 40 as compared to 25 GeV 2 ). 
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Figure Caption 

Fig. 1 The 7p — > T(lS)p cross section as a function of the 7^ centre-of-mass energy W. The 
data points are the ZEUS |TJ and preliminary HI measurements. The horizontal lines 
indicate the range of W sampled by the measurements, the vertical lines the systematic 
and statistical errors added linearly. The continuous curves are obtained from the two 
QCD calculations (model I and II) described in the text. The dotted curve corresponds 
to the naive leading-order formula (0), the dashed curve shows our prediction for method 
II if we neglect the corrections from the real part of the amplitudes. 
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